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A control axis referenced model of the NASA/JPL 70-m antenna structure is com- 
bined with the dynamic equations of the servo components to produce a comprehensive 
state variable (matrix) model of the coupled system. An interactive Fortran program for 
generating the linear system model and computing its salient parameters is described. Re- 
sults are produced in a state variable , block diagram , and in factored transfer function 
forms to facilitate design and analysis by classical as well as modern control methods. 


I. Introduction 

The upgrade of the NASA/JPL 64-m antennas to 70-m 
apertures has added considerable mass and inertia moment on 
top of the existing alidade with resultant decreases in the natu- 
ral frequencies of the structure. Because the combined compli- 
ances of the alidade and gear reducers separate the autocolli- 
mator from the antenna servo motors and tachometers, the 
increased inertia can affect the dynamics of autocollimator 
based pointing control. A study was undertaken to assess 
the impact on the axis servos of the increased inertia and 
decreased frequencies and to provide a more complete model 
for the new servo design. This article describes the methodol- 
ogy of combining the dynamics of the structure with those 
of the servo and cites results for both the 64-m and 70-m 
antennas. 

The development of a condensed antenna structural model 
[1] provides a control axis referenced representation of the 
structure dynamics in a compact form. When integrated with 
the dynamic properties of the hydraulic actuators and the con- 
trol electronics, the result is a more comprehensive model for 


design, analysis, and simulation of the axis servos. The com- 
posite model is derived by coupling the linear differential 
equations describing the structure with those of the control 
components to produce a state variable model. 

The structure model consists of a relatively large residual 
inertia to which are coupled individual modal inertias. For the 
elevation axis, two modes of the alidade structure are added. 
The elasticity of the control actuators, the compliance and 
inertia of the hydraulic system, the drive motors, and the 
characteristics of the servo compensation networks are also 
superimposed. In the most inclusive of the optional forms 
available, the elevation model representing five modal inertias 
and two alidade modes results in a nineteenth order linear sys- 
tem model. 

The antenna pointing system employs two modes of posi- 
tion feedback derived from either a shaft angle encoder driven 
through a precision gear reducer from the bullgear, or an opti- 
cal autocollimator mounted at the rear of the apex of the main 
reflector. The autocollimator mirror is attached to an hour 
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angle-declination mount on a pedestal isolated from that of 
the antenna. The structure model includes the transformation 
coefficients relating the displacement at the autocollimator 
reference structure to the arbitrary coordinates of the indi- 
vidual inertias. These coefficients enable system modeling of 
displacements at the autocollimator as well as at the encoder 
and tachometers, all of which are connected by nonzero 
flexibilities. 

The state variable representation has one limitation in that 
system response properties are not evident by inspection when 
large system matrices are involved. For example, the existence 
of a pole in the right-half s plane, evident in a root locus dis- 
play, is not recognizable by inspection of a linear system 
matrix unless the matrix happens to be in a form where the 
eigenvalues are recognizable. The matrix form of system repre- 
sentation does, however, provide great flexibility to accommo- 
date various computer processing methods. Thus, both trans- 
fer function and block diagram representations are readily 
derived from the open loop as well as the closed rate loop 
system matrices. The results are available in a variety of forms 
suitable for classical as well as modern control system design 
and analysis. 

An interactive Fortran program was developed for generat- 
ing the model in state variable (matrix) form. Operating on an 
IBM PC, the program provides options for adjustment of the 
complexity of the model and of various parameter values. It 
displays the resultant linear system matrix and computes the 
corresponding system poles from the eigenvalues of the system 
matrix. Zeros of each of the tachometer, encoder, and autocol- 
limator responses are determined from the eigenvalues of par- 
titioned matrices formed according to the Mason gain formula 
[2] . Zeros of the autocollimator response are computed from 
the weighted sum of encoder and individual mode zeros where 
the weighting factors are the transformation coefficients men- 
tioned earlier. This indirect method was used as an expedient 
to avoid coding an algorithm for the more general solution of 
zeros of the transformed matrix. 

II. Structure Dynamics Model 

The form of the condensed structure mode for the eleva- 
tion axis is shown schematically in Fig. 1. In Fig. 
and 0 B correspond, respectively, to the angular rates of 
motion at the gear reducer attachment, the reducer output 
pinion (normalized with respect to the pinion to bullgear 
ratio), and the twin bullgears. The combined stiffness of the 
four gear reducers (two driving each bullgear) is represented by 
a single spring, K G , shown connected between the pinion and 
the “single” bullgear. The angular rates of the individual modal 
inertias 6 } through 0 N have an indirect correspondence to 
the motion of the physical antenna. A set of dimensionless 


coefficients, alpha, relate the motion of the Intermediate 
Reference Structure to through d N . 

While the azimuth axis of the antenna employs a single, sta- 
tionary bullgear and moving gear reducers, its dynamic motion 
can be described using the same model as in the elevation axis. 
A small error results, however, because in azimuth, the reducer 
housing rotates with the alidade, such that in 360 degrees of 
azimuth motion each azimuth drive motor shaft makes N-l 
revolutions relative to the pedestal and only N revolutions 
relative to the alidade. The error arises because the motor 
torque-speed characteristic is referenced to alidade coordinates, 
while acceleration torque is proportional to motion in inertial 
coordinates. However, since the gear ratio, N , is large, the 
resulting error is negligible. 

Using the coordinate definitions of Fig. 1, the equations 
of motion for the individual inertias,^., are derived. In general, 
a small damping, £>-, is associated with each spring, K { . 

J.6. + D.(d. - 0 B ) + K.{d. -S B ) = 0 for / = 1 to N 

( 1 ) 

j J b *( k c*T,Wb-'Z, m-v*w = 

1 1 


J A2 e A2+<. K Al +K A2>°A2- K AX 9 Al = 0 W 

III. Hydraulic Motor Model 

A schematic representation of the hydraulic system and 
gear reducers is shown in Fig. 2. Each gear reducer includes 
two fixed displacement, axial piston hydraulic motors driving 
a spur gear reducer. One motor provides control and the other 
produces a constant value countertorque. The countertorque 
motors are coupled to a regulated hydraulic pressure supply 
and connected so as to produce opposite torque outputs from 
two of the four reducers. This arrangement eliminates back- 
lash in all gear meshes except the single mesh between each 
control motor and the first intermediate gear. The counter- 
torque has the added benefit of increasing the incremental 
torque stiffness of the reducer geartrains. Control motor 
modulation is provided by a four-port hydraulic servovalve 
interconnecting a regulated high pressure hydraulic source and 
return with the motors. 
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The dynamic model of the drive motors and hydraulic sys- 
tem is derived according to established practices. Because the 
high gain of the rate servo loop diminishes the overall effects 
of friction, leakage, and valve pressure variations, a compara- 
| tively simple, second order, linearized model can be justified. 

Neglecting leakage, the oil flow, Q m , through the hydraulic 
^ motor is given by 

| Q = VO 

! m 

where V is the volumetric displacement of the motor and 6 m is 
| the rotational speed at the output shaft. 

Equating the fluid input power and the mechanical output 
power yields the relationship between motor torque, T m , and 
hydraulic pressure , P 


T = PV 

m 

The oil flow through the hydraulic servo valve is described 
by the equation of flow through a sharp edge orifice of vari- 
able size 


3,=Wv 1/2 (S) 

where P v is the pressure drop across the valve, is the flow 
constant of the valve, and I v is the valve coil current. 

I 

The valve pressure drop equals the pressure difference 
between the regulated supply and the motor. Since both analy- 
tic and field test results show the steady state pressure drop is 
| relatively constant in the range of 2000 to 2500 psi (13790 to 
17238 kPa), a piecewise linear approximation of the above 
orifice equation is sufficient. Transient pressure drops caused 
| by large acceleration or high wind torque will effectively 
decrease the rate loop gain and increase valve damping relative 
to the values represented by the linear approximation. With 
the aid of a root locus diagram, it can be shown that the net 
effect on the rate loop stability is negligible. 

| The linearized valve equation combined with, the motor and 
compressible flow equations is illustrated in block diagram 
form in Fig. 3 where C and J m represent the hydraulic compli- 
ance and motor inertia, respectively. The valve gain and damp- 
ing constants and/), respectively, are derived by partial dif- 
ferentiation of the valve flow Eq. (5) with values of P v and Q v 
that represent mean operating conditions. The damping param- 
eter D may be increased to include motor leakage and other 
equivalent sources of damping. 

In rearranging the above equations to a form compatible 
with the structure equations, a provision for the motion of 
the flexible structure supporting the reducer is required. This 


is accomplished by equating the motor rotation to the differ- 
ence between 6 m and 0 Al in Fig. 1. Thus by application of the 
Mason gain rule to the diagram of Fig. 3 and with the defini- 
tion of the reducer natural frequency 


CO 



1/2 


the gear reducer equation is obtained as: 
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' a+ £ S + "«)<*»-*4> + 


NJ 

r m 




2 
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( 6 ) 


where Q = K p I v is the effective no-load valve flow and N r is 
the overall motor-to-axis gear ratio. 


The large value of gear ratio N r (28700) permits the omis- 
sion of the acceleration torque term that arises from the differ- 
ence between reducer and inertial coordinates, with a small 
error resulting. 


IV. System Equations 

Equating the reducer shaft torque, T x , to the torque trans- 
mitted by the reducer stiffness, K G completes the equations 
defining the coupled actuation-structure system as follows: 

-V < 7 > 

To accommodate analyses of the servo rate loops as well as 
the position loops, and also to facilitate computation of trans- 
fer function parameters, the system equations are developed 
in two steps. First, the open loop linear system matrix and 
associated input and output vectors are derived from Eqs. (1) 
through (7) above. Subsequently, the closed loop system 
matrix is formed by augmentation of the open loop matrix by 
inclusion of the rate feedback and compensation gains. The 
system state equations are derived from Eqs. (1) through (7). 
above using the state variable definitions in Table 1. The vari- 
able definitions are generalized to accommodate either the azi- 
muth or elevation axis and a variable number of structure 
modes, N. Because the azimuth structure model is based upon 
the assumption of a rigid pedestal, state variables correspond- 
ing to $ A1 and 2 are absent from the azimuth model. 

It will be seen that with the exception of the third, (2iV+5)th, 
and (2Af+9)th variables, all rates and accelerations are relative 
to stationary coordinates. The third variable corresponds to 
the effective rotor-stator rotation of the hydraulic motor and 
of the tachometer. The (2/V+5) and (2N+9) variables corre- 
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spond to the effective extension of the gear reducer torsional 
stiffness for the azimuth/elevation axes, respectively. This tor- 
sional extension is related to the reducer output torque 
through the stiffness parameter K G . 

The resulting state equations for the elevation axis are listed 
in generalized form in Table 2. Using 


[F]X+ [G] U 

( 8 a) 

= [H)X 

( 8 b) 


with X the state vector, X the time derivative of X, U the sys- 
tem input, and Y the output, the corresponding linear system 
matrices are listed in Table 3. The system block diagram is 
shown in Fig. 4. The azimuth equations, matrices, and block 
diagram are similar to their elevation counterparts, except the 
four alidade states are omitted. 


V. Derivation of Transfer Functions 

The linear transfer functions relating the inputs and outputs 
are expressed as ratios of factored polynomials in the Laplace 
operator, s. The respective numerator and denominator factors 
are of the form (5 - zero) and (5 - pole) and are related to the 
system matrices through the traditional equations expressing 
the conditions for the poles and zeros of the response. 

For system poles (denominator factors): 


corresponding matrices representing the equations of the 
graphs. Using this equivalence, the transfer function numerator 
thus becomes the determinant of the partitioned matrix 
formed by deleting from the system equations those state 
variables included along the forward path in a flow graph 
representation of the equations. Since forward paths are recog- 
nizable in system equations, this partitioning can be accom- 
plished without actually constructing the graph. The num- 
erator factors are thus determined from the eigenvalues of the 
partitioned matrix formed above. 

The zeros for the autocollimator and elevation encoder out- 
put responses are computed by superposition of individual 
components of the respective responses. This superposition 
avoids both a coordinate transformation of the system matrix 
and also the relatively complicated application of the Mason 
rule to a transformed matrix. Because it involves weighted 
summations of characteristic polynomials with subsequent 
factoring, this method is subject to numerical inaccuracy as 
matrix size increases. Good accuracy has resulted for models 
including three structure modes. The accuracy of this method 
could be improved and its usefulness extended to more modes 
by frequency scaling in such a way as to reduce the numeric 
range of the polynomials. Using superposition, the general 
case for the encoder response characteristic polynomial,/^, is 
given by 


K r K 

p = — p p 

E J B B W A 


(10b) 


[sJ-F] = [0] (9) 

For system zeros (numerator factors): 


1 

^3 

1 


~X(s) 

H 0 


_U(s)_ 


where the response zeros are the roots of P E \P B and P A are 
the characteristic polynomials of the bullgear and alidade 
responses, respectively. They are the characteristic polynomials 
of submatrices formed according to the method described 
above. Equation (10b) can be applied to the azimuth axis by 
equating P A to zero. 


Equation (9) is an eigenvalue equation and the system poles 
are the eigenvalues of the linear system matrix. Software for 
eigenvalue evaluation is available for DOS microcomputers. 
For the more complicated case of Eq. (10) for evaluating the 
zeros, a more convenient method based on the Mason gain rule 
was employed. An advantage of this alternate method is that 
it reduces the dimensions of the matrices to be processed, thus 
improving numerical accuracy. 

The method is based upon an adaptation of Mason’s [2] 
signal flow graph gain (transfer function) determination, to 
state variable representations. Mason relates transfer function 
denominator and numerator to “determinants” of flow graphs 
and of certain subgraphs. It can be shown that these determi- 
nants of flow graphs are identical to the determinants of 


The complexity of the elevation autocollimator response is 
reduced by partitioning into two smaller systems with the 
results combined so as to avoid processing polynomials of high 
order. The validity of this simplification is evident from the 
bullgear response submatrix formed by deleting the state vari- 
ables included in the forward path between input and the bull- 
gear rate. This submatrix is formed from the F matrix of 
Table 3 by deletion of rows and columns 1 through 4. Since 
this submatrix consists of two diagonal blocks, with one block 
representing the alidade and the other representing the tipping 
structure, the eigenvalues are simply the combination of those 
of the individual blocks. The complete solution of the autocol- 
limator response is the weighted superposition of the responses 
of the bullgear and of the individual inertias, all of which are 
derived from this block diagonal matrix. The zeros of the 
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autocollimator response are thus the roots of characteristic 
polynomial P AC defined by Eqs. (10c) and (lOd): 


(12) 



matrices having the same subscripts, is the bullgear response 
submatrix described above, and the F)- matrices are derived 
from F bb by row, column deletion by the aforementioned 
principle. Coefficients a 0 • ■ • a N in Eq. (10c) are the trans- 
formation coefficients provided in the structure model. 


The physical significance of this partitioning is explained by 
considering the individual uncoupled responses of the alidade 
and the tipping structure to a torque input applied at the ele- 
vation gear reducer. The frequencies of infinite compliance 
(i.e., zero stiffness) of the alidade are frequencies of bullgear 
response zeros and are not influenced by the tipping structure. 
Additional frequencies of zero bullgear response are the natu- 
ral frequencies of the individual tipping structure modes and 
are not influenced by the presence of the alidade. The parti- 
tioning is thus valid and no loss of generality results from its 
use. 


VI. Closed Rate Loop Model 

Modeling of the closed rate loop configuration is accom- 
plished by extending the open loop equations and matrices to 
include the rate feedback and loop compensation. In both the 
64-m and 70-m rate loops, tachometer based feedback, a lag / 
lead network, and a lead/lag network are employed to obtain a 
high degree of stiffness at low frequency and a comparatively 
narrow noise bandwidth. Since the 50 Hz servovalve band- 
width is sufficiently wide to have negligible effect on loop 
dynamics, the rate feedback can be modeled by the two phase 
compensation networks and a gain parameter. 

The transfer function of the tachometer-network-amplifier, 
servovalve cascade is: 

Q<s ) A :,r fr + Z,)(i + Z,) 

The two real poles of Eq. (11) result in two additional 
states with the following additional state equations: 


X M ~ K r ( Z 1 P l ^ x 3 P 1 X M 
X M+ 1 = 2 ~ P l) X 3 + ( Z 2 ~ P 2 ~ P 2 X M+ 1 

(13) 

where 

M = 2 A+6 for azimuth and 

M = 27V+10 for elevation 

For the closed loop case the equation forx 4 becomes 

*4 = - + ~f X A +X m +X m + l ( 14 > 

The closed rate loop linear system matrix is obtained by 
augmenting the open loop matrix of Table 3 according to 
Eqs. (12), (13), and (14) and is shown in Table 4. 

The various transfer functions for the open and closed loop 
systems are compiled according to the following properties of 
linear systems: 

(1) all transfer functions in a given system have identical 
poles. 

(2) the zeros in any given transfer function are invari- 
ant with respect to gain changes in other parts of 
the system not touching the forward path under 
consideration. 

Property (1) implies that the tachometer, encoder, and 
autocollimator transfer functions all have identical poles, thus 
eliminating a need for repeated computations. Because the 
loop closing gain is in a feedback path, property (2) allows 
the use of zeros computed for the open loop case in the closed 
loop transfer functions. The real zero at -2.2, appearing in the 
closed loop functions, results from an electronic compensation 
network in the forward path. 

VII. Numerical Results 

A generalized Fortran program for generating the linear sys- 
tem matrix and computing the corresponding poles and zeros 
was developed and executed on an IBM PC. The program pro- 
vides options for selection of axis, display and adjustment of 
parameter values, and adjustment of the number of structure 
modes included. As the number of modes is adjusted down- 
ward, the inertia of the rejected modes is added to the residual 
inertia, thus maintaining accuracy of the total inertia. Results 
are written to disk-files in a format compatible with postpro- 
cessing programs. 
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The Fortran program was used to derive models for the azi- 
muth and elevation axes of both 64-m and 70-m configurations. 
Hardware parameter values available for the four configura- 
tions are listed in Table 5. In some cases the number of struc- 
ture modes modeled was reduced from the number available 
in order to reduce computation and data space. For the 64-m 
azimuth case the 4th mode was eliminated, for 70-m azimuth 
the 5th mode was eliminated, and for 70-m elevation the 3rd, 
4th, and 5th modes were merged into a weighted composite. 
The best estimate of structure damping ratio available is 0.003, 
which was observed in factory tests of the 70-m quadripod. 
Because of computation error in some of the zeros computa- 
tions resulting from nonzero damping, all computations were 
run with zero structure damping. The introduction of appro- 
priate damping would displace the complex zeros from the 
imaginary axis and cause a similar displacement of the com- 
plex poles. However, since the poles are already damped by 
the rate loop, the introduction of structure damping would 
have small effect on the overall results. 

The poles and zeros of the open and closed rate loop trans- 
fer functions are listed in Table 6. In cases where inertia values 
other than those of Table 5 were used, the actual values are 
included in Table 6. Frequency response plots of amplitude vs. 
radian frequency for the 70-m azimuth and elevation axes are 
shown in Figs. 5 and 6. 

The differences between the tachometer, encoder, and 
autocollimator responses are due to the flexibility of the gear 
reducers and structure between the respective devices. In the 
elevation autocollimator responses, the low frequency (2 Hz 


for 64-m, 1.5 Hz for 70-m) resonant peak and subsequent roll- 
off result from the flexibility of the alidade. It occurs because 
alidade deflections in elevation are sensed by the autocollima- 
tor but not by the encoder. Both the encoder and autocol- 
limator responses for both 64-m axes show a high frequency 
roll-off beginning at 6.7 to 7.2 Hz. For the 70-m configuration 
the roll-off frequencies are nearly identical, 6.5 to 6.8 Hz. This 
result is explained by the fact that most of the inertia increase 
resulting from the upgrade is associated with the first three 
modes while the residual inertias are relatively unchanged. 

VIII. Conclusions 

The models described have enabled the design of the rate 
and position servos for the 70-m configuration with a mini- 
mum of on-site adjustment. The azimuth axis structure model, 
see Table 5, includes three low frequency modes at 1.27 to 
1.42 Hz, a frequency roughly half that of the lowest mode in 
the 64-m configuration. These modes are a cause for concern 
because of their low frequencies, close spacing, and relatively 
large associated inertias which make damping by the servo 
loop difficult to achieve. The presence of these modes necessi- 
tated an increase of the motor shaft mounted inertia wheels as 
a means of diminishing their effect on the control system. 

The results summarized in Table 6 provide improved defini- 
tion of the encoder and autocollimator response characteristics 
as compared with those based on the assumption of a “rigid” 
alidade structure. These results will be useful in future efforts 
to improve autocollimator based pointing. They do not neces- 
sitate immediate changes as a result of the upgrade. 
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Table 2. State equations for elevation 


Table 1. State variable definitions 


Variable 

Symbol 

Description 


e B 

Bullgear angular rate 

*2 

*1 

Bullgear acceleration 

*3 

9 M~ 9 Al 

Motor/tach rate 

X 4 


Hydraulic torque/ 

X 2i+3 

6. 

t 

Angular rate at the inertia J. 
for / = 1 to N 

X 2i+4 

X 2i+3 


X 2/V+S 

9 m~ 6 b 

Motor-bullgear angle difference, 
azimuth only 

X 2N+ 5 

e A 1 

Angular rate at alidade 1, elevation 
only 

X 2N+6 

X 2N+S 

Elevation only 

X 2N+1 

9 A 2 

Angular rate at alidade 2, elevation 
only 

X 2N+8 

X 2N+1 

Elevation only 

X 2N+9 

9 M~ 9 B 

Motor-bullgear angle difference, 
elevation only 


*1 = X 2 


IK „ JL K. 


N 


'.‘St? 


« T •b/ 


AT. K g K g 

7- x 2i + 3 + 7Z x 3 + ~7T x 2N + s 

l D D 


X 3 *4 J X 2N+9 
M 


7 D , Q 2 

X. = -w ■ x,, +*=r co 

4 m 3 £ 4 F m 


for i = 1 to TV 


• . 2 „ °i 2 °i 

X 2i+4 w / X i * J. X 2 X 2/+3 J. *2/+4 


K c k g *g +k ai k ai 

X 2N + 6- J a1 X 1- J a1 X 3- J m X 2N+S+^ X 2N+2 


A1 


K A1 + K A2 


X 2N+S T~ X 2N+S 77; X 2N+1 

A2 .42 


X 2N+9 X l +X 3 + X 2N+S 


N = number of structure modes in model 
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Table 3. Linear system matrices for open loop elevation axis 
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Table 3. (contd) 


OUTPUT VECTORS, H 

FOR TACHOMETER OUTPUT RESPONSE 


H t =[0 0 1 0 0] 

FOR ENCODER OUTPUT RESPONSE 

H Ea = I 1 0 0 0 0 1 Azimuth only 

H Ee = [ 1 0 0 0 ... -1 .... 0] Elevation only 


FOR AUTOCOLLIMATOR MOUNT RESPONSE 

H ac = [ aO 0 0 0 al 0 a2 0 a3 0 a4 0 a5 0 . . . . 0 ] 

NOTES: 

All input and output vectors have lengths equal to the dimension of the corresponding linear system matrix, F. 

The -1 in the Elevation Encoder response is in column 2N + 5, where N is the number of structure modes modeled. 

Coefficients al. . .a5 are defined as a part of the structure model. When fewer than five structure modes are modeled the a’s correspond- 
ing to unmodeled modes are replaced with zero and aO is increased an equal amount such that the sum of aO + al + . . . aN = 1 where 
N is the number of structure modes modeled. 


Table 4. Linear system matrices for closed rate loop 



where T is a matrix with 1.0 in Row 4, Col. 3 and zeros elsewhere 



Hj, H£, are formed from their open loop counterparts by 
adding two zero elements following the last element 
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Table 5. 64- and 70-m antenna structure and servo parameter values 



64-m AZ 

70-m AZ 

64-m EL 

70-m EL 

RESIDUAL (Base) INERTIA, J B 

0.1783 

0.1813 

0.0840 

0.1066 

(0.2418) 

(0.2459) 

(0.1139) 

(0.1446) 

INERTIA Ratio, J X /J B 

0.3097 

0.2356 

0.3782 

1.1746 

INERTIA Ratio, ^/J B 

0.1092 

0.4368 

0.3936 

0.2383 

INERTIA Ratio, J 3 /J B 

0.0505 

0.2494 

0.2729 

0.0433 

INERTIA Ratio, J 4 /J fi 

0.0975 

0.0436 


0.0812 

INERTIA Ratio, J 5 /J B 


0.0132 


0.0121 

MOTOR INERTIA, ^ 

0.6640 

1.0000 

0.6640 

0.6640 

(0.9006) 

(1.3563) 

(0.9006) 

(0.9006) 

ALIDADE INERTIA, J A1 



0.0145 

0.0145 



(0.0197) 

(0.0197) 

ALIDADE INERTIA, 



0.1330 

0.1330 



(0.1804) 

(0.1804) 

FREQ of (Kgear/Jg) 0 ' 5 

38.3639 

38.0399 

66.5865 

59.1247 

FREQ, mode 1, R/s 

15.0500 

7.9670 

19.3830 

14.7910 

FREQ, mode 2, R/s 

25.7600 

8.3120 

20.7600 

17.6930 

FREQ, mode 3, R/s 

32.0800 

9.8900 

26.4020 

20.2760 

FREQ, mode 4, R/s 

36.0200 

13.7470 


21.8340 

FREQ, mode 5, R/s 


16.9020 


25.6290 

FREQ of Hyd Motor 

9.5940 

7.8000 

9.5940 

9.5940 

RE part, Motor freq, D/C 

1.2000 

1.2000 

1.2000 

1.2000 

FREQ of Alidade 1 



67.8258 

67.8258 

FREQ of Alidade 2 



32.0237 

32.0237 

IRS Transformation coeff. al 

0.1318 

0,1331 

0.1421 

0.2939 

a2 

0.2022 

0.2767 

0.1722 

0.0977 

a3 

0.1040 

0.1169 

0.2090 

0.0266 

a4 

0.0618 

0.0099 


0.0309 

a 5 


0.0376 


0.0050 

RATE LOOP GAIN, K R /to 2 
R m 

2302 

2302 

2302 

2302 


RATE LOOP COMPENSATION NETWORKS, all axes 


(s + 2.2) (s + 7.1) 
(s + 0.12) ’ (s + 81) 


NOTES: 

64-m and 70-m azimuth models are for 90 degree elevation. 

Inertia values are ft-lb-s 2 , (kg-M 2 ) and are referred to motor shaft. 
Natural frequencies are in radians/s. 

Coefficient aO = 1.0000 - al - a2 - a3 - a4 - a5. 


Table 6a. Transfer functions for 64-m azimuth axis: 3 structure modes modeled 


(FILE: AZ6415) 

RESIDUAL (Base) INERTIA, JB 0.1956 (0.2653) 


INERTIA Ratio, Jl/JB 
INERTIA Ratio, J2/JB 
INERTIA Ratio, J3/JB 
INERTIA Ratio, J4/JB 
INERTIA Ratio, J5/JB 
MOTOR INERTIA, JM 

OPEN RATE LOOP 


VALVE CURRENT TO MOTOR SHAFT TRANSFER FUNCTION 
POLES ZEROS 


Real 

Imag 

Real 

Imag 

-5.0933E-Q03 

4.4378E+001 

O.OOOOE+OOO 

4.0634E+001 

-5.0933E-003 

-4.4378E+001 

0.0000E+000 

-4.0634E+001 

-1.7687E-003 

3.1629E+001 

O.OOOOE+OOO 

3.1012E+001 

-1.7687E-003 

-3.1629E+001 

O.OOOOE+OOO 

-3.1012E+001 

-4.2904E-003 

2.5545E+001 

O.OOOOE+OOO 

2.4712E+001 

-4.2904E-003 

-2.5545E+001 

O.OOOOE+OOO 

-2.4712E+001 

-2.3554E-002 

1.5505E+001 

O.OOOOE+OOO 

1.4625E+001 

-2.3554E-002 

-1.5505E+001 

O.OOOOE+OOO 

-1.4625E+001 

-5.6529E-001 

7.8390E+000 



-5.6529E-001 

-7.8390E+000 




CLOSED RATE LOOP 

RATE LOOP INPUT TO TACHOMETER 
RATE TRANSFER FUNCTION 


POLES ZEROS 


Real 

Imag 

Real 

Imag 

-3.1927E+000 

4.4833E+001 

O.OOOOE+OOO 

4.0634E+001 

-3.1927E+000 

-4.4833E+001 

O.OOOOE+OOO 

-4.0634E+001 

-5.5539E-001 

3.1564E+001 

O.OOOOE+OOO 

3.1012E+001 

-5.5539E-001 

-3.1564E+001 

O.OOOOE+OOO 

-3.1012E+001 

-8.6159E-001 

2.5261E+001 

O.OOOOE+OOO 

2.4712E+001 

-8.6159E-001 

-2.526 1E+001 

O.OOOOE+OOO 

-2.4712E+001 

-5.3122E-001 

1.4465E+001 

O.OOOOE+OOO 

1.4625E+001 

-5.3122E-001 

-1.4465E+001 

O.OOOOE+OOO 

-1.4625E+001 

-1.2952E+0G1 

1.6448E+001 



-1.2952E+G01 

-1.6448E+001 



-1.4483E+G00 

O.OOOOE+OOO 

-2.2000E+000 

O.OOOOE+OOO 

-4.4686E+001 

O.OOOOE+OOO 

-8.1000E+001 

O.OOOOE+OOO 


0.2821 

0.0995 

0.0460 

0.0000 

0.0000 

0.6640 (0.9006) 

CLOSED RATE LOOP 

RATE LOOP INPUT TO ENCODER RATE TRANSFER FUNCTION 


POLES ZEROS 


Real 

Imag 

Real 

Imag 

-3.1927E+000 

4.4833E+001 



-3.1927E+000 

-4.4833E+001 



-5.5539E-001 

3.1564E+001 

O.OOOOE+OOO 

3.2080E+001 

-5.5539E-001 

-3.1564E+001 

O.OOOOE+OOO 

-3.2080E+001 

-8.6159E-001 

2.5261E+001 

O.OOOOE+OOO 

2.5760E+001 

-8.6159E-001 

-2.5261E+001 

O.OOOOE+OOO 

-2.5760E+001 

-5.3122E-001 

1.4465E+001 

O.OOOOE+OOO 

1.5050E+001 

-5.3122E-001 

-1.4465E+001 

O.OOOOE+OOO 

-1.5050E+001 

-1.2952E+001 

1.6448E+001 



-1.2952E+001 

-1.6448E+001 



-1.4483E+000 

O.OOOOE+OOO 

-2.2000E+000 

O.OOOOE+OOO 

-4.4686E+001 

O.OOOOE+OOO 

-8.1000E+001 

O.OOOOE+OOO 


CLOSED RATE LOOP 


RATE LOOP INPUT TO AUTOCOLLIMATOR 
RATE TRANSFER FUNCTION 

POLES ZEROS 


Real 

Imag 

Real 

Imag 

-3.1927E+000 

4.4833E+001 



-3.1927E+000 

-4.4833E+001 



-5.5539E-001 

3.1564E+001 

O.OOOOE+OOO 

3.6678E+001 

-5.5539E-001 

-3.1564E+001 

O.OOOOE+OOO 

-3.6678E+001 

-8.6159E-001 

2.5261E+001 

O.OOOOE+OOO 

2.8297E+001 

-8.6159E-001 

-2.5261E+001 

O.OOOOE+OOO 

-2.8297E+001 

-5.3122E-001 

1.4465E+001 

O.OOOOE+OOO 

1.5985E+001 

-5.3122E-001 

-1.4465E+001 

O.OOOOE+OOO 

-1.5985E+001 

-1.2952E+001 

1.6448E+001 



-1.2952E+001 

-1.6448E+001 



-1.4483E+000 

O.OOOOE+OOO 

-2.2000E+000 

O.OOOOE+OOO 

-4.4686E+001 

O.OOOOE+OOO 

-8.1000E+001 

O.OOOOE+OOO 
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Table 6b. Transfer functions for 70-m azimuth axis: 4 structure modes modeled 


(FILE: AZ7028) 


RESIDUAL (Base) INERTIA, JB 

0.1837 

INERTIA Ratio, Jl/JB 

0.2325 

INERTIA Ratio, J2/JB 

0.4311 

INERTIA Ratio, J3/JB 

0.2462 

INERTIA Ratio, J4/JB 

0.0430 

INERTIA Ratio, J5/JB 

0.0000 

MOTOR INERTIA, JM 

1.0000 


OPEN RATE LOOP CLOSED RATE LOOP 


VALVE CURRENT TO MOTOR SHAFT TRANSFER FUNCTION 

RATE LOOP INPUT TO ENCODER RATE TRANSFER FUNCTION 

POLES 

ZEROS 

POLES 

ZEROS 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

-3.1307E-0Q3 

4.2053E+001 

O.OOOOE+OOO 

3.8858E+001 

-1.5156E+000 

4.2566E+001 



-3.1307E-003 

-4.2053E+001 

O.OOOOE+OOO 

-3.8858E+001 

-1.5156E+000 

-4.2566E+001 



-2.8774E-003 

1.3785E+001 

O.OOOOE+OOO 

1.3706E+001 

-4.8021E-002 

1.3696E+001 

O.OOOOE+OOO 

1.3747E+001 

-2.8774E-003 

-1.3785E+001 

O.OOOOE+OOO 

-1.3706E+001 

-4.8021E-002 

-1.3696E+001 

O.OOOOE+OOO 

-1.3747E+001 

-6.1942E-002 

1.0306E+001 

O.OOOOE+OOO 

9.8101E+000 

-6.7226E-002 

9.7472E+000 

O.OOOOE+OOO 

9.8900E+000 

-6.1942E-002 

-1.0306E+001 

O.OOOOE+OOO 

-9.8101E+000 

-6.7226E-002 

-9.7472E+000 

O.OOOOE+OOO 

-9.8900E+000 

-7.6929E-002 

8.9059E+000 

0.00G0E+000 

8.2325E+000 

-3.6465E-002 

8.1685E+000 

O.OOOOE+OOO 

8.3120E+000 

-7.6929E-002 

-8.9059E+000 

O.OOOOE+OOO 

-8.2325E+000 

-3.6465E-002 

-8.1685E+000 

O.OOOOE+OOO 

-8.3120E+000 

-5.9236E-003 

8.063 1E+000 

O.OOOOE+OOO 

7.9107E+000 

-9.1012E-002 

7.8285E+000 

O.OOOOE+OOO 

7.9670E+000 

-5.9236E-O03 

-8.063 1E+000 

O.OOOOE+OOO 

-7.9107E+000 

-9.1012E-002 

-7.8285E+000 

O.OOOOE+OOO 

-7.9670E+000 

-4.491 IE-001 

6.1701E+000 



-8.401 1E+000 

1.3145E+001 



-4.491 IE-001 

-6.1701E+000 



-8.401 1E+000 

-1.3145E+001 







-6.0554E+001 

O.OOOOE+OOO 

-8.1000E+001 

O.OOOOE+OOO 





-1.4472E+000 

O.OOOOE+OOO 

-2.2000E+000 

O.OOOOE+OOO 


CLOSED RATE LOOP 



CLOSED RATE LOOP 


RATE LOOP INPUT TO TACHOMETER 

RATE LOOP INPUT TO AUTOCOLLIMATOR 


RATE TRANSFER FUNCTION 



RATE TRANSFER FUNCTION 


POLES 

ZEROS 

POLES 

ZEROS 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

-1.5156E+000 

4.2566E+001 

O.OOOOE+OOO 

3.8858E+001 

-1.5156E+000 

4.2566E+001 



-1.5156E+000 

-4.2566E+001 

O.OOOOE+OOO 

-3.8858E+001 

-1.5156E+000 

-4.2566E+001 



-4.8021E-002 

1.3696E+001 

O.OOOOE+OOO 

1.3706E+001 

-4.802 IE-002 

1.3696E+001 

O.OOOOE+OOO 

1.3353E+001 

-4.8021 E-002 

-1.3696E+001 

O.OOOOE+OOO 

-1.3706E+001 

-4.8021E-002 

-1.3696E+001 

O.OOOOE+OOO 

-1.3353E+001 

-6.7226E-002 

9.7472E+000 

O.OOOOE+OOO 

9.8101E+000 

6.7226E-002 

9.7472E+000 

O.OOOOE+OOO 

1.3009E+001 

-6.7226E-002 

-9.7472E+000 

O.OOOOE+OOO 

-9.8101E+000 

-6.7226E-002 

-9.7472E+000 

O.OOOOE+OOO 

-1.3009E+001 

-3.6465E-002 

8.1685E+000 

O.OOOOE+OOO 

8.2325E+000 

-3.6465E-002 

8.1685E+000 

O.OOOOE+OOO 

9.2957E+000 

-3.6465E-002 

-8.1685E+000 

O.OOOOE+OOO 

-8.2325E+000 

-3.6465E-002 

-8.1685E+000 

O.OOOOE+OOO 

-9.2957E+000 

-9.1012E-002 

7.8285E+000 

O.OOOOE+OOO 

7.9107E+000 

-9.1012E-002 

7.8285E+000 

O.OOOOE+OOO 

8.0628E+000 

-9.1012E-002 

-7.8285E+000 

O.OOOOE+OOO 

-7.9107E+000 

-9.1012E-002 

-7.8285E+000 

O.OOOOE+OOO 

-8.0628E+000 

-8.401 1E+000 

1.3 145E+001 



-8.4011E+000 

1.3145E+001 



-8.401 1E+000 

-1.3145E+001 



-8.401 1E+000 

-1.3145E+001 



-6.0554E+001 

O.OOOOE+OOO 

-8.1000E+001 

O.OOOOE+OOO 

-6.0554E+001 

O.OOOOE+OOO 

-8.1000E+001 

O.OOOOE+OOO 

-1.4472E+000 

O.OOOOE+OOO 

-2.2000E+000 

O.OOOOE+OOO 

-1.4472E+000 

O.OOOOE+OOO 

-2.2000E+000 

O.OOOOE+OOO 
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Table 6c. Transfer functions for 64-m elevation axis: 3 tipping structure and 2 alidade modes modeled 


(FILE: EL6414) 


RESIDUAL (Base) INERTIA, JB 

0.0840 (0.1139) 

INERTIA Ratio, Jl/JB 

0.3782 

INERTIA Ratio, J2/JB 

0.3936 

INERTIA Ratio, J3/JB 

0.2729 

INERTIA Ratio, J4/JB 

0.0000 

INERTIA Ratio, J5/JB 

0,0000 

MOTOR INERTIA, JM 

0.6640 (0.9006) 


OPEN RATE LOOP 


CLOSED RATE LOOP 


VALVE CURRENT TO MOTOR SHAFT TRANSFER FUNCTION RATE LOOP INPUT TO ENCODER RATE TRANSFER FUNCTION 
POLES ZEROS POLES ZEROS 

Real Imag Real Imag Real Imag Real Imag 


-2.2606E-005 
-2.2606E-005 
5.2814E-004 
5.2814E-004 
3.8028E-005 
3.8028E-005 
2.9061E-004 
2.9061E-004 
1 .7055E-004 
1.7055E-004 
9.0286E-G02 
9.0286E-002 
5.0871E-001 
5.0871E-001 


1.8661E+002 

-1.8661E+002 

4.3688E+001 

-4.3688E+001 

3.1688E+001 

-3.1688E+001 

2.3835E+001 

-2.3835E+001 

1.9984E+001 

-1.9984E+001 

1.4827E+001 

-1.4827E+001 

8.0624E+000 

-8.0624E+000 


O.OOOOE+OOO 

0.0000E+000 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 


1.8528E+002 
-1.8528E+002 
4.3309E+001 
-4.3309E+001 
3.1678E+001 
-3. 1678E+001 
2.3806E+001 
-2.3806E+00I 
1.9976E+001 
-1.9976E+001 
1.2709E+001 
-1.2709E+001 


-3.2473E-002 
-3.2473E-002 
-2.9969E-001 
-2.9969E-001 
-9.6810E-003 
-9.6810E-003 
— 2.2261E-002 
-2.226 IE-002 
-4.3576E-003 
-4.3576E-003 
-4.9794E-001 
-4.9794E-001 
-2.0885E+001 
-2.0885E+001 
-3.7370E+001 
-1.4486E+000 


1.8668E+002 

-1.8668E+002 

4.3696E+001 

-4.3696E+001 

3.1684E+001 

-3.1684E+001 

2.3809E+001 

-2.3809E+001 

1.9975E+001 

-1.9975E+001 

1.2338E+001 

-1.2338E+001 

1.6166E+001 

-1.6166E+00I 

O.OOOOE+OOO 

O.OOOOE+OOO 


O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 


-8.1000E+001 

-2.2000E+000 


4.4681E+001 

-4.4681E+001 

3.1582E+001 

-3.1582E+001 

2.3968E+001 

-2.3968E+001 

1.9799E+001 

-1.9799E+001 

1.3135E+001 

-1.3135E+001 


O.OOOOE+OOO 

O.OOOOE+OOO 


CLOSED RATE LOOP 


CLOSED RATE LOOP 


RATE LOOP INPUT TO TACHOMETER 
RATE TRANSFER FUNCTION 


RATE LOOP INPUT TO AUTOCOLLIMATOR 
RATE TRANSFER FUNCTION 


POLES 


ZEROS 


POLES 


ZEROS 


Real Imag Real Imag 


Real Imag Real Imag 


-3.2473E-002 
-3.2473E-002 
-2.9969E-001 
-2.9969E-001 
-9.6810E-003 
-9.6810E-003 
-2.2261E-002 
-2.2261E-002 
.3576E-003 
.3576E-003 
.9794E-001 
.9794E-001 
-2.0885E+001 
-2.0885E+001 
-3.7370E+001 
-1 .4486E+000 


1.8668E+002 
-1.8668E+002 
4.3696E+001 
-4.3696E+001 
3.1684E+001 
-3.1684E+001 
2.3809E+001 
-2.3809E+001 
1.9975E+001 
-1.9975E+001 
1.2338E+001 
-1.2338E+001 
1.6166E+001 
—1.6 166E+00 1 
O.OOOOE+OOO 
O.OOOOE+OOO 


O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 


-8.100GE+001 

-2.2000E+000 


1.8528E+002 
-1.8528E+002 
4.3309E+001 
-4.3309E+001 
3.1678E+001 
-3.1678E+001 
2.3806E+001 
-2.3806E+001 
1.9976E+001 
-1 9976E+001 
1.2709E+001 
-1.2709E+001 


O.OOOOE+OOO 

O.OOOOE+OOO 


-3.2473E-002 

-3.2473E-002 

-2.9969E-001 

-2.9969E-001 

-9.6810E-003 

-9.6810E-003 

-2.2261E-002 

-2.2261E-002 

-4.3576E-003 

-4.3576E-003 

-4.9794E-001 

-4.9794E-001 

-2.0885E+001 

-2.0885E+001 

-3.7370E+001 

-I.4486E+000 


1.8668E+002 

-1.8668E+002 

4.3696E+001 

-4.3696E+001 

3.1684E+001 

-3.1684E+001 

2.3809E+001 

-2.3809E+001 

1.9975E+001 

-1.9975E+001 

1.2338E+001 

-1.2338E+001 

I.6166E+001 

-1.6166E+001 

O.OOOOE+OOO 

O.OOOOE+OOO 


O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 

O.OOOOE+OOO 


-8.1000E+001 

-2.2000E+000 


7.2267E+001 

-7.2267E+001 

3.4188E+001 

-3.4188E+001 

3.0055E+001 

-3.0055E+001 

2.2681E+001 

-2.2681E+001 

1.9844E+001 

-1.9844E+001 

O.OOOOE+OOO 

O.OOOOE+OOO 
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Table 6d. Transfer functions for 70-m elevation axis: 3 tipping structure and 2 alidade modes modeled 


INERTIA AND FREQ OF 3rd MODE ARE ADJUSTED COMPOSITES OF MODES 3,4,5 


(FILE: EL7012) 


RESIDUAL (Base) INERTIA, JB 

0.1066 (0.1446) 

INERTIA Ratio, Jl/JB 

1.1746 

INERTIA Ratio, J2/JB 

0.2383 

INERTIA Ratio, J3/JB 

0.1362 

FREQ, mode 3, R/s 

21.7400 

INERTIA Ratio, J4/JB 

0.0000 

INERTIA Ratio, J5/JB 

0.0000 

MOTOR INERTIA, JM 

0.6640 (0.9006) 


OPEN RATE LOOP CLOSED RATE LOOP 


VALVE CURRENT TO MOTOR SHAFT TRANSFER FUNCTION 

RATE LOOP INPUT TO ENCODER RATE TRANSFER FUNCTION 

POLES 

ZEROS 

POLES 

ZEROS 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

-2.3597E-005 

1.8433E+002 

O.OOOOE+OOO 

1.8299E+002 

-3.3846E-002 

1.8441E+002 



-2.3597E-005 

-1.8433E+002 

O.OOOOE+OOO 

-1.8299E+002 

-3.3846E-002 

-1.8441E+002 



-6.7820E-004 

4.1970E+G01 

O.OOOOE+OOO 

4.1540E+001 

-3.5693E-001 

4.1960E+001 

O.OOOOE+OOO 

4.2889E+001 

-6.7820E-004 

-4.1970E+G01 

O.OOOOE+OOO 

-4.1540E+001 

-3.5693E-001 

-4.1960E+001 

O.OOOOE+OOO 

-4.2889E+001 

-1.3184E-004 

2.7651E+001 

O.OOOOE+OOO 

2.7629E+001 

-2.0608E-002 

2.7636E+001 

O.OOOOE+OOO 

2.7644E+001 

-1.3184E-004 

-2.765 1E+001 

O.OOOOE+OOO 

-2.7629E+001 

-2.0608E-002 

-2.7636E+001 

O.OOOOE+OOO 

-2.7644E+001 

-7.1847E-004 

2.0753E+001 

O.OOOOE+OOO 

2.0711E+001 

-2.5872E-002 

2.0708E+001 

O.OOOOE+OOO 

2.0737E+001 

-7.1847E-004 

-2.0753E+001 

O.OOOOE+OOO 

-2.071 1E+001 

-2.5872E-002 

-20.708E+001 

O.OOOOE+OOO 

-2.0737E+001 

-1.8263E-003 

1.7085E+001 

O.OOOOE+OOO 

1.7039E+001 

-1.5653E-002 

1.7033E+001 

O.OOOOE+OOO 

1.7050E+001 

-1.8263E-003 

-1.7085E+001 

O.OOOOE+OOO 

-1.7039E+001 

-1.5653E-002 

-1.7033E+001 

O.OOOOE+OOO 

-1.7050E+001 

-1.8590E-001 

1.2952E+001 

O.OOOOE+OOO 

9.8583E+000 

-3.2123E-001 

9.4835E+000 

O.OOOOE+OOO 

1.0185E+001 

-1.8509E-001 

-1.2952E+001 

O.OOOOE+OOO 

-9.8583E+000 

-3.2123E-001 

-9.4835E+000 

O.OOOOE+OOO 

1.0185E+001 

-4.1066E-001 

7.1233E+000 



-2.0812E+001 

1.6437E+001 



-4.1066E-001 
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Fig. 1. Structure model schematic representation 



Fig. 2. Countertorque system, 64/70-m antenna 



Fig. 3. Hydraulic motor linearized block diagram 




















